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Abstract 

The 3+1 GRMHD equations for Schwarzschild spacetime in Rindler coordinates with isothermal state of 
plasma are formulated. We consider the cases of non-rotating and rotating backgrounds with non-magnetized 
and magnetized plasmas. For these cases, the perturbed form of these equations are linearized and Fourier 
analyzed by introducing plane wave type solutions. The determinant of these equations in each case leads 
to two dispersion relations which give value of the wave number k. Using the wave number, we obtain 
information like phase and group velocities etc. which help to discuss the nature of the waves and their 
characteristics. These provide interesting information about the black hole magnetosphere near the horizon. 
There are cases of normal and anomalous dispersion. We find a case of normal dispersion of waves when the 
plasma admits the properties of Veselago medium. Our results agree with those of Mackay et al. according 
to which rotation of a black hole is required for negative phase velocity propagation. 

Keywords : 3+1 formalism, GRMHD, Rindler spacetime, isothermal plasma. 

1 Introduction 

General Relativity (GR) is a beautiful scheme for describing the gravitational field and the corresponding 
equations. This theory is believed to apply to all forms of interaction whether these are between the charged 
particles or large scale gravitational structures, including compact objects like black holes. Black holes are still 
mysterious 1 . While physicists have been grappling the theory of black holes, astronomers have been searching 
for real- life examples of black holes in the universe [2] . It is proved that black holes exist on the basis of study of 
effects which they exert on their surroundings. With their enormous gravitational fields, they greatly affect the 
surrounding plasma medium (which is highly magnetized) and hence plasma physics in the vicinity of a black 
hole has become a subject of interest in astrophysics. In the immediate vicinity of a black hole, GR applies and 
it is therefore of interest to formulate plasma physics problem in this context. 

In large scale gravitational structures, one often encounters magnetic fields. All the massive stars, neutron 
stars and black holes carry energy flux which produces a relatively large magnetic field. In order to understand 
such phenomena, we have a developed theory of magnetized plasma called theory of magnetohydrodynamics 
(MHD). In this theory, perfect MHD flow is important which is the flow of plasma with negligible charge sepa- 
ration and with a magnetic field frozen into it [5^. Since all the compact objects have strong gravitational fields 
near their surfaces 4 , it is important to study the general relativistic effects on physical processes taking place 
in their vicinity. The MHD theory with the effects of gravity is called general relativistic magnetohydrodynam- 
ics (GRMHD). The GRMHD equations include Maxwell's equations. Ohm's law, mass, momentum and energy 
conservation equations. These help us to study stationary configurations and dynamic evolution of conducting 
fluid in a magnetosphere. They are also required to investigate various aspects of the interaction of relativistic 
gravity with plasma's magnetic fleld. 

The general relativistic physics of black holes and plasmas takes place in 3+1 dimensional space while the 
corresponding relativistic laws are written in three-dimensional language. This formalism is well-suited to carry 
non-relativistic intuition of physicists about plasmas, hydrodynamics and stellar dynamics into the arena of 
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black holes and GR. The formalism (also called ADM formalism) was originally developed by Arnowitt et 
al. [5] and was motivated by several startling results proved in 1970 using the black hole viewpoint 
Thorne and Macdonald [T2]-[I1] extended the formulation to electromagnetic fields of the black hole theory. 
The wave propagation theory in the Friedmann universe was investigated using 3+1 formalism by Holcomb and 
Tajima [T3], Holcomb [TO] and Dettmann et al. [TT]. Khanna 181 derived the MHD equations describing the 
two component plasma theory of the Kerr black hole in this split. Anton et al. [19] used the formalism to 
investigate various test simulations and discussed magneto-rotational instability of accretion disks. Komissarov 
PP] discussed Blendford-Znajec monopolc solution using 3+1 formalism in black hole electrodynamics. Zhang 
[21]- [22] formulated the black hole theory for stationary symmetric GRMHD using this formalism with its 
applications in Kerr geometry. 

The gravitational perturbations away from the Schwarzschild background have been discussed by Regge and 
Wheeler [23] . This treatment was extended by Zerilli 24] who showed that the perturbations corresponding to 
a change in mass, the angular momentum and charge of the Schwarzschild black hole are well-behaved. The 
decay of non well-behaved perturbations has been investigated by Price 25J. The quasi-static electric problem 
was solved by Hanni and Ruffini [8J who proved that the lines of force diverge at the horizon for the observer 
at infinity. Wald (26j derived the solution for electromagnetic field occurring when a stationary axisymmetric 
black hole is placed in a uniform magnetic field aligned along the symmetry axis of black hole. A linearized 
treatment of plasma waves for special relativistic formulation of the Schwarzschild black hole was developed by 
Sakai and Kawata [27]. Using 3+1 ADM formalism, Buzzi et al. [28 extended this treatment to waves in general 
relativistic two component plasma propagating in radial direction. They investigated the one dimensional radial 
propagation of transverse and longitudinal waves close to the Schwarzschild horizon. In a recent paper, Sharif 
and Umber [29] have found some interesting properties of cold plasma in the vicinity of a Schwarzschild black 
hole horizon using 3+1 GRMHD equations. They have also verified that the rotation is a cause of negative 
phase velocity propagation which coincides with the result of Mackay et al. [30 . 

The aim of this paper is to discuss the properties of waves in isothermal plasma. We develop and expound 
a theoretical method relevant to the study of relativistic fluid model used in applications to astrophysics and 
plasma physics. We consider Schwarzschild spacetime in Rindler coordinates which gives an example of the 
essential ideas of the horizon and the 3+1 split without distracting complications of curved spacial three ge- 
ometry. This investigation would lead to support the viewpoint that information from a black hole cannot be 
extracted. 

The GRMHD equations in 3+1 formalism using isothermal plasma in the Schwarzschild magnetosphere are 
obtained. The component form of the equations is extracted by involving linear perturbations in each of the 
cases of either non-magnetized or magnetized, non-rotating or rotating perfect fluid. The condition for the plane 
harmonic wave is imposed on all the perturbations. The equations are then Fourier analyzed and solved to get 
the dispersion relations with the help of software Mathematica. Two dispersion relations are found in each case 
and the wave numbers are deduced uniquely in terms of angular frequency. The wave velocity, group velocity, 
refractive index and the corresponding quantities give the properties of plasma near the horizon expressed by 
the perfect fluid. These properties are shown in terms of graphs. 

The scheme of the paper is as follows. The next section contains basic informatory results for 3+1 formalism. 
The results for the GRMHD equation are given and energy equation is obtained. These equations are restricted 
for the case of perfect fluid with perfect MHD flow assumption for Rindler approximation of the Schwarzschild 
spacetime. The equations arc further modified by adding the condition of isothermal plasma. Section 3 is 
devoted to the case of non-rotating background for which we find the dispersion relations. Section 4 is furnished 
with the values of k for rotating non-magnetized plasma. In section 5, we generalize the problem to the case of 
rotating magnetized plasma. In the last section, we shall summarize and discuss the results. 

2 3+1 Spacetime Modeling 

The astrophysical theory in early 1960s was dominated by the frozen star viewpoint [3T]. As long as the 
viewpoint prevailed, physicists failed to realize that the black hole can be dynamical, evolving, energy-releasing 
object. Many problems of current interest in black-hole astrophysics cannot be posed in frozen-star viewpoint, 
because the viewpoint has difficulty producing unambiguous boundary conditions at the event horizon. Such 
problems where it is inadequate to replace the black hole horizon, in picture and calculations, by a surface of 
no return, a new viewpoint for the black hole physics was developed - the membrane paradigm (a marriage 
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of horizon formalism with 3+1 formahsm). This is mathematically equivalent to the standard, full, general 
relativistic theory of black holes for the physics outside the horizon. Although this formalism is not treated 
in standard relativity textbooks, it has been developed in great detail by many researchers and has played an 
important role 14j. 

In 3 + 1 formalism, the line clement of the spacetime can be written as [22] 

ds^ = -a^dt^ + jij{dx' + I3^dt){dx^ + /3^dt), 

where a is the lapse function, /3* are the components of shift vector and 7^ are the components of the spatial 
metric. All these quantities are functions of the coordinates t, . A natural observer, associated with this 
spacetime called fiducial observer (FIDO), has four- velocity n perpendicular to the hypersurfaces of constant 
time t. Notice that we are using geometrized units. 

The Rindler approximation of the Schwarzschild line element is [M] 

ds^ ^ -a^{z)dt^ + dx^ + dy^ + dz^ (2.1) 

which is an analogue of the Schwarzschild spacetime with z, x and y as radial r, axial (j) and poloidal 6 directions 
respectively. In Rindler coordinates, the lapse function a is "^here ru is the radius of the Schwarzschild 
black hole. This function clearly vanishes at the horizon which we can place at z = 0. It is useful to think of 
this absolute space with fiducial observer which never moves from its fixed location. The relationship between 
the universal time t and FIDO measured time r can be expressed by the time lapse a = The Schwarzschild 
black hole is non-rotating, hence the shift vector f3 — 0. This is the minimum configuration of the Kerr black 
hole out of which no energy can be extracted. 

2.1 Perfect GRMHD Equations 

In magnetohydrodynamic treatment, the plasma is represented as a perfect fluid. The perfect MHD flow 
assumption is given by 

E + VxB = 0, (2.2) 

where V, B and E are the velocity, magnetic and electric fields of the fluid measured by FIDO. This assumption 
shows that fluid possesses no electric fleld in rest-frame and the electric fleld is perpendicular to the magnetic 
fleld. 

It is known that unlike cold plasma approximation, isothermal plasma contains plasma pressure which cannot 
be ignored. The pressure gradient term is to be included. In this model, the energy equation is also introduced 
but the term involving the heat flux is neglected. Because of the vanishing thermal conductivity, the plasma is 
non- viscous and consequently the non-diagonal terms of the pressure dyad are all equal to zero. FIDO measured 
local energy conservation law |21j is given by 

^+0e + ^Wi£a'n - - Av.(a2s) + iy/? : W + E.j. (2.3) 

Here e, S, W^^ and j represent the mass energy density, energy flux, stress tensor and current vector of the 
electromagnetic fleld in three-dimensions respectively, 9 is the expansion rate of the FIDO's four- velocity, £t is 
the time derivative along shifting congruence (Lie derivative with respect to global time in a usual fashion) and 
^ = ^(^ — /3-V) is the rate of change of a three-dimensional vector which lies in the absolute space according 
to the FIDO. When we substitute the following values for perfect fluid [TJj 

e - (a*Po-p(1-V2))72, 
S = fipaj^V, 
W = /ipo7^V (g) V + 

with fi the speciflc enthalpy, pa is the rest- mass density, p is the pressure, and perfect MHD condition, Eq. (|2.3p 
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takes the following form 

Po7'^ + A^7'^ + 2poAi7^V.^ - ^ + (a^pot' - pW + 2poM7'V.a + poM7'(V.V) 

W7'V.(V.V)/3 +p(V./3)} + -L{poA^7'^'^^' + P7^-''}i^t7., + 7^[(V x B).(V x (aB)) 

+a(V X B).-^(V X B) + (V X B).(V x B.V)/3 + ea{V x B).(V x B)] = 0, (2.4) 
dr 

where 

For the Schwarzschild black hole, FIDO measured Faraday's law, equation of evolution of the magnetic field, 
local conservation laws of mass and momentum are given by Eqs.(2.17)-(2.20) ^9]. These equations along with 
local conservation law of energy Eq. (|2.4p become 

— = V X (aV X B), (2.5) 

— + (aV.V)B - (B.V)(aV) + (V.(aV))B = 0, (2.6) 
ot 

+ {(aV).V}(poM) + PoM7'V.^ + poM7'V.(aV.V)V + pofi{V.{aV)} = 0, (2.7) 

+P^l'V. (-^ + (V.V)m| = -P0M7'«. -P,. + ^(V X B),V.(V X B) - ^(aB)^ 
{a ot J 47r 87ra^ 

+ -^(aS,).,B^' - -^[B X {V X (V X (aV x B))}]„ (2.8) 



-poP7'(V.V) + 

47ra 



(V X B).(V X (aB)) + (V X B). — (V x B) 

ot 



= 0, (2.9) 



where a is the gravitational acceleration. Notice that here Eq. (|2.4p is the new equation in addition to Eqs.(2.17)- 
(2.20) [21]. These constitute the perfect GRMHD equations for fluid with non-zero pressure in the vicinity of 
a black hole horizon. 



2.2 GRMHD Equations for Isothermal Plasma in Schwarzschild Spacetime 

The isothermal equation of state can be expressed by the following equation [21] 

fi = ^ ^ ^ — constant. (2-10) 

PO 

Here p is the mass density and pressure of the fluid. 

For the Rindler approximation of Schwarzschild geometry, FIDOs four-velocity expansion rate and shift 
vector vanish. Thus the perfect GRMHD equations in 3+1 formalism for isothermal state of plasma turn out 
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to be 



— =Vx(aVxB), (2.11) 

— + (aV.V)B - (B.V)(aV) + BV.(aV) = 0, (2.12) 
^^^tZl + (aV.V)(p + p) + {p + Ph^^-^ + {p + ph^V.{aV.V)V + (p +p)V.(aV) = 0, (2.13) 

iP + Ph' + IJ) % + ip + ph'v.v, ^s.s,} (^1 + v.v) 

- ^B.S,) V^kV' = -(P + P)7'«. -E. + ^V.(V X B)(V X B), 
(aB)2 {aB,)^jBi 1 



87ra^ 47ra 47ra 



[B X {V X (V X (aV X B))}]„ (2.14) 



7'(-4 + V.V)(p + p) + 2(p + p)7'V.(l| + V.V)V + 2ip + p)j^V.a - 
a ot a ot a at 

+ {p + p)-1^V.V + -^[(V X B).(V X (aB)) + (V x B).|-(V x B)] = 0. (2.15) 

Ana at 

Eas. (|2.1ip - (|2.15p are the perfect GRMHD equations for isothermal plasma in the vicinity of the Schwarzschild 
black hole. In the rest of the paper, the Fourier analysis method will be applied to the perturbed form of these 
equations (including the restrictions of backgrounds) to obtain the dispersion relations. 

3 Non- Rotating Background 

In non-rotating background, we apply the Fourier analysis method to the perturbed GRMHD equations given 
in Appendix A. These equations will lead to the dispersion relations. We obtain solution of dispersion relations 
as graphs. The same procedure will be applied to rotating background. 

We assume the plane wave type solutions of the form g-^i^t-kz) perturbations (Appendix A). Thus the 
perturbed variables can be expressed as follows 

p{t, z) = cie-'('^*-'=^), p{t, z) = C2e-'("*-'=-), v,(t, z) = Cje-'^^^-'^-), 6,(<, z) = C4e-'('^*-'=^), (3.1) 

where k is the wave number and oj is the angular frequency. 

Using Eq. (|3.ip the Fourier analyzed form of the perturbed Eqs. (jA7P - (|All[) specified for the non-rotating 
background, it follows that 

C4 = 0, 3.2 

a 

ikci = 0, (3.3) 
ci{— ptcj + ikpau — (uap)' — au^^'^pu'} + C2{—p^^ + tkpau + (uap)' + au^j'^pu'} 

+C3{p + p) |a(l + -i'^u^)ik - a(l - 2-i'^u'^){l + j'^u'^)— - icuj'^u^ = 0, (3.4) 
„2r„ I „,„,'n I „,2„,2m I „ r„„,2r„ , „,„,//i , „,2„,2\ 



CiP7 {az + uu'{l +'-f u )} + C2{p7 {a^ -I- uu'{l + 7"u )} + ikp + p'} 



+C3{P + P)l^ 



a 



(1 + -I'^u^) + mfc + {u'{l + 7^u^)(l + 47^^^) + 27^0 J 



P7^ V iku] + Yup' + p2Yua^ + p(l + 27^u^)7^u' \ ci 



= 0, (3.5) 



— p{l^ ^ 1) + i-k^f^up + "f^up + 2'f'^upaz + p(l + 27^u^)7^u' > C2 
a J 

+ {p + ph^ I ^^7^w + (ifc + il^uu' + a^)(l + 2j^u^) - — I C3 = 0. (3.6) 
{a u ) 

Eqs. (|3.2p and (|3.3p yield that C4 = 0, i.e., the magnetic field has no effect of black hole gravity as well as time. 
We would like to mention here that Eqs. (|3.4p - p.6p also give non-magnetized plasma. 
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Numerical Solutions 

We consider that the magnetosphere is filled with stiff ffuid for which p — constant — p. We also assume that 
the time lapse a = z. Using these values, the mass conservation law in 3-dimensions, i.e., a{p+p)^u — constant 
gives u = ^^^^'i ■ Due to the magnetic field parallel to the wave number and the use of the longitudinal part of 
momentum equation, the isothermal plasma shows only the longitudinal (electron or ion) plasma waves. 

When we use these values, the determinant of the coefficients of constants ci, C2 and C3 in Eqs. p.4p - 
(j3.6p leads to a complex dispersion equation. On comparing the real and imaginary parts we obtain two 
dispersion equations. The real part shows an equation of the type Ai{z)k'^ + A2{z,uj)k + ^3(2,0;) = which 
on solving gives two values of the wave number k. The imaginary part gives the dispersion equation of the 
type Ai{z)k^ + ^2(2, Lo)k^ + A3{z,u)k + ^4(2, w) = which provides three (including two complex conjugate) 
values of the wave number. Using the values of k, the quantities like phase velocity {vp = j), refractive index 
(n = ^), its change with respect to angular frequency (^) and group velocity [vg = = ^^^dn ) can be 

obtained. These quantities help us to investigate the properties of the waves as well as their dispersion. 

The two dispersion relations obtained from the real part give the same value of k. The corresponding graphs 
are given in Figure 1. The wave number obtained from the imaginary part is shown in Figure 2. 

In Figure 1, the wave number k becomes infinite at 2: = which means that the waves vanish at the horizon 
due to effect of immense gravity. No wave with zero angular frequency exist which is a physical consequence. 
The wave number decreases as we go away from the horizon which means that the waves lose energy as their 
distance from the horizon increases. The wave number increases with the increase in angular frequency. The 
phase and group velocities are positive and increase when we go away from the horizon. We can deduce that the 
increase in gravity decreases the wave velocity as well as the speed with which energy travels. Both the phase 
and group velocities are equal. The refractive index is greater than one and there is no change with respect to 
angular frequency hence the dispersion is not normal there [32]. 

Figure 2 shows that the wave number is infinite at z = 0. In the region 0<z<l, 0<a;<6x 10^^, the 
wave number takes some complex value and the waves are evanescent there, otherwise the wave number k is 
positive. It decreases as we move away from the event horizon of the Schwarzschild black hole. The increase 
in angular frequency increases the wave number. The phase and group velocities are equal except in the region 
< w < 6 X 10^^. Since the phase velocity increases as we go away from the horizon, the waves gain energy 
as they go far away from the event horizon of the black hole. Refractive index is positive and greater than one 
but its derivative with respect to angular frequency has random variation. Hence the medium is not of normal 
dispersion. 

4 Rotating Non-Magnetized Background 

The perturbed form of GRMHD Eqs. (|2.1ip - (|2.15p meant for rotating non-magnetized background are given in 
Appendix B. For the sake of Fourier analysis, we assume that the perturbations (Appendix B) are of the plane 
wave form, i.e., e-«('^*-'=^) . Thus the perturbed variables take the form 

p{t, z) = cie-^("*-'=-), p{t, z) = c2e-'('^*-'=^\ v^, z) = cae-'^'^*"'^^), v.,{t, z) - C4e-^(-*-'=-). (4.1) 

Using Ea. (|4.ip . we can write the Fourier analyzed Eqs. (|B7P - (jB10l) for rotating non- magnetized background 
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Surface Plot 



of Wave Number 





Figure 1: The waves decrease as they move away from the horizon. The phase and group velocities are equal. 
The dispersion is not normal. 
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Surface Plot 



of Wave Number 




Figure 2: The wave number decreases as the waves move away from the event horizon. In most of the region 
the phase and group velocities are equal and the dispersion is not normal. 
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as follows 



Ci{{—lll) + bakujp — a'up — au'p — aup — a^'^up{VV' + uu')} 
+C2{(— + Laku)p + a'up + au'p + aup' + aj'^up{VV' + uu')} 



-Luj-f^u + Lkail + j^u^) - a <^ (1 - 2-/^u^)(l + -f^u'^) 2j\^VV' 

' u 



+C4{p + p)[-/^V{-LLu + ikau) + a-f'^u{{l + 2-f^V'^)V' + 2-^'^uVu'}] = 0, 
cip-i^u{{l + 7^2)^' + -i^uVu'} + C2P7^ii{(l + 7^2)1/' + -f^uVu'} 

-bUJ 



(4.2) 



c?,{p + p)l 



-i^uV 



+ t/sM + {(1 + 272y2)(i + 272u2) „ -yV^}!/' + 272uVu'(l + 272u2) 



+C4(p + p)7^ 



(1 + 7^2) + Lk-^ + 7^u{(l + ^-i^u^)uu + 4(1 + 7^2)^^'} 



= 0, 



(4.3) 



c\pT{a^ + um'(1 + 7^?r) + 7^u^yF'} + C2\pT{az + W(l +7^u^) + 7^u^T/T/'} + p' + Lkp\ 



C3{p + Ph^ 

+C4(p + p)7'' 



(1 + 7^^^) + t/cuj + u'{l + j^u^)il + 47V) + 2m72{(1 + 2-f^u^)VV' + aj 

uV ( — + Lku] + mV(1 + 47^^^) + 2y{(l + 2j^u^)uu' + aj = 0, 



(4.4) 







' —LLO 




Cl 















+ tfcu ) + 2p7^u{az + 7^(FV^' + uu')} + p7^u' + p'7^" 



u 



C2 



-(7^ - 1) + ifc7^M 



+2p7^u{a^ +7^(\/y' + uu')} +P7^u' +p'7^u] +cz{p + p){ 2j'^u \-j^Lk{l + 2'fu^) 

y a 

-7^— + ^-f'u^iyV + uu') + -i^iyV + uu') + 27''uu' + 7^0^(1 + 2-i'^u^)\ 
u J 



+C4(p + p) 



27V 



+ t/cu + 27''{372uF(V^F' + uu') + uV + uT/a^)} 



= 0. 



(4.5) 



Numerical Solutions 



We consider time lapse a = z and assume that = u in the stiff fluid of constant density. Using mass 



conservation law in three dimensions (with constant mass density), we obtain the value of u 



Since 



the magnetic field is absent, the plasma contains the longitudinal electron plasma waves, ion plasma waves or 
the transverse electromagnetic waves. 

The determinant of the coefficients of constants ci, C2, C3 and C4 from Eqs. (|4.2p - (|4.5p gives a complex 
dispersion equation. The real part leads to dispersion equation of the type ^1(2;) A;'* + ^2(2, u))kr' + ^3(2, LS)k^ + 
A^{z^uj)k + A'^(z^lS) = 0. From the imaginary part, we obtain equation of the type A\(z)k^ + A2{z,u!)k'^ + 
A3{z,Lu)k + A4{z,uj) = 0. The dispersion relation from the real part gives four values of k out of which two 
are interesting. The other two values turn out to be imaginary in the whole region. The imaginary part of the 
determinant gives three values out of which two are complex conjugates. The wave numbers obtained from the 
real part are shown in Figures 3 and 4 while Figure 5 demonstrates the value of k obtained from the imaginary 
part. 

We see from Figure 3 that fc is positive for < z < 10, 0.011 < < 10 while it is negative for < z < 
1, < < 0.011. It is interesting to note that in this small region, the phase velocity and refractive index 
are negative whereas the group velocity is positive. The fluid shows the properties of Veselago medium |33| . 
The wave number decreases with the increase in z and increases with the increase in angular frequency in the 
region 1 < 2; < 10, 0.011 < < 10. In the neighborhood of w = 0, the phase velocity is large. It decreases with 
the increase in z with the exception of the interval where it has complex values. The group velocity shows the 
similar behavior with the exception that it is complex at w = 0. In the region 0<z<2, 2 < lo < 10, the index 
of refraction is greater than one and also ^ > 0. The dispersion is therefore normal in this small region but 
not normal otherwise. 

In Figure 4, the wave number becomes imaginary in the region < w < 0.5 and hence evanescent waves 
exist there. The wave number increases with an increase in lo and the waves lose energy as we go away from 
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Surface Plot 



of Wave Number 




Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 3: The wave number decreases when waves move away from the event horizon. A small region is showing 
properties of Veselago medium. The dispersion is not normal in most of the region. 
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Surface Plot of Wave Number 




Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 4: The wave number decreases as the waves move away from the event horizon. The dispersion of waves 
is found to be not normal. 
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the horizon. In the region < z < 10, 0.5 < a; < 10, the wave number as well as the phase and group velocities 
are real. The refractive index is greater near the horizon due to immense gravitational field and decreases away 
from the horizon. The change of refractive index with respect to the angular velocity is negative in the whole 
region. Thus the dispersion is not normal. 

Figure 5 gives the description of the wave number obtained from the imaginary part of the determinant. 
The wave number increases as they move towards the horizon of the black hole. It also increases as the angular 
frequency of waves increases. The waves are evanescent in the region < z < 0.112 while the real waves lie in 
the region 0.112 < z < 10. The phase velocity is zero at z — 0.112 and it increases towards 1 as z increases. The 
refractive index is greater than one in this region but its change with respect to angular frequency is negative 
which shows that dispersion is not normal. 

5 Rotating Magnetized Background 

This is the most general case. The perturbed GRMHD equations for this situation are given in Appendix C. Since 
these equations include the additional magnetic field, therefore the Fourier analysis procedure demands that the 
magnetic field perturbations should have harmonic space and time dependence, i.e., of the form b ^ ^-i-i^^t-kz) ^ 
These perturbed variables admit the following notations 

6,(i,z) = C5e-'("*-'=^), 6,(i,z) = C6e-'(-*-'=^). (5.1) 

When we substitute these values in component form of perturbed Eas. (jC2[) and (|C3[) (Appendix C), it follows 
that 

C5 y h ikuj = 0, t/cc5 — 0. 

Both these equations yield that C5 = and hence = 0. Substituting this value in addition to the values from 
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Surface Plot 



of Wave Number 




Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 5: The wave number decreases as the waves move away from the horizon of the black hole. The waves 
are not normally dispersed. 
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Eqs. (|4.ip and (|5.ip in Eqs. ljCip . (|C4p - (|C7p . we obtain the following Fourier analyzed equations 



— csKaA)' + ikaX} + 04(0:' + ika) — ceKau)' — luj + tkau} = 0, 
ci {p {~LUJ + ikau) — a'up — au'p — aup' — aupj'^ {VV' + itw')} 
+C2 {p (— iti' + ikau) + a' up + au'p + aup' + aupj'^ {VV' + uu')} 



-Luj-f^u + Lka{l + j'^u^) -a\{l- 2j^u^){l + -f^u'^) 2-/'^u^VV' 



+ C3{p + p) 



+Ci{p + p)j^ [{-LLu + Lkau)V + au{{l + 2j^V^)V' + 2-/^uVu'}] = 0, 
cipj^u{{l + "/^V^)V' + -/^uVu'} + C2p-/^u{{l + -/^V^)V' + -/^uVu'} 



+C3 



An 



B^u. 



+ iP + Ph^ {{(1 + 2-f^u^)il + 27^2) - -f'^V^}V' + 272(1 + 2-f^u^)uVu'} + ^(Aa)' 



-C4 



{p+p)r{i+Yv') + 



47r J a 

+ ip + p)7^u{(l + Aj'^V^)uu' + 4(1 + -f^V^)VV'} 
B^ 



+ Lku{ {p+p)"f^{l+"f'^V^) 

B'^ua'' 



B^ 

An 



Ana 



An 



cq < {1 — u'^)Lk H (1 — u^) — uu 



0, 



ci7V[a^ + u{{l+^'^u^)u' +-i'^VuV'}] + C2[-i'^p{a:, + m{(1 + 72^^)^' + 7^^-^^'}} + ikp + p'] 



'C3 



- Up + Ph''{i + i^u^) + 



An 



+ < (p + p)7 (1+7 w )- 



An 



iku 



+ \{p + P)7'{"'(1 + 7'«')(1 + 47'u2) + 2u72{(l + 272^^)^^' + a,}} - ^^(aA)' 

47ra 



+C4 



- I (p + p)7^^F _ ^ I ^ + I (p + p)^^uV +^Yk'. 

+ + p)7^{uV(l + 47^2) + 2y{a^ + W(l + 27V)}} + 

(Aa)' 



AB^a'u 

Ana 



B 

H C6 { A(l - u2)tA: + A — (1 - u^) - Ami' + 



An 
ci7^ 



0, 



p + ifcu^ + 2pu{a^ +72(FV^' + um')} + pu' + p'u + C2 p | -^(7^ - 1) + 6/172^ 



+2p7"u{a2 + 7 {VV' + uu)} + p^ u + p"f tt] + C3 



(p + p) <^ 27*u h 7^ifc(l + 27^u^) 



-7^— + 676^2(^1/' + uu') + 7^(1/^' + uw') + 27'*um' + 7^0^(1 + 2-i'^u^) 
u 



2i\p + p){V 



a 



iku 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



B 

H {A(aA)' - ■u(q;A)'(uA -V) - L{kau + uj)\{u\ -V) + C4 

Ana 

+{-i-i'^uV{VV' + uu') + uV + uFa^)}) + {-(aA)' + a'uiuX ~ F) + (t/cau + tw)(uA - V)} 

Ana 

B^ 

+C6- \u{a\)' + [a' - u{au)' + tfc(l - u^)}{uX - V)] = 0. 

Ana 

Numerical Solutions 

We take the same assumptions as in the case of rotating non-magnetized plasma, i.e., u — 
tuting this value in = i + Au of 21 , it turns out that A = 1 — 



vW2 



(5.6) 



= V. Substi- 



This shows that the magnetic field 
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diverges near the horizon (as discussed in J3 ). Also, is taken to be 

The matrix of the coefficients of constants ci, C2, C3, C4 and cg is of order 5x5 whose determinant gives a 
dispersion equation quartic in k from the real part and a quintic from the imaginary part. The real part leads 
to four interesting solutions whereas the imaginary part is difficult to solve to get exact solutions. It is obvious 
that the magnetic field is non-zero and the wave number is in arbitrary direction to the magnetic field. Since 
the external magnetic field is present, the plasma is anisotropic. This magnetic field has profound effects in 
plasma wave modes discussed in the previous section. 

The real part gives dispersion relation of the type Ai{z)k^+A2{z, uj)k^+A3{z, uj)k'^+A4{z, uj)k+A5{z, uj) — 0. 
It gives four values of k out of which two arc complex conjugate roots and are not interesting. The two interesting 
roots provided by the real part are shown in Figures 6 and 7. 

In Figure 6, waves are evanescent in the regions < z < 1.65 and < < 2.5 x 10"'*. The wave number 
becomes infinite at zero angular frequency as well as at z = and hence waves do not exist there. Real waves 
are present in rest of the region. We see that the wave number is large near the horizon and decreases gradually 
afterwards as z increases. Also, k increases with an increase in the angular frequency. The phase and group 
velocities increase as we go away from the horizon. Also, the phase velocity is zero at zero angular frequency 
but the group velocity is infinite there. The refractive index is less than one in the whole region and hence the 
dispersion is not normal. 

Figure 7 shows that the evanescent waves lie in the regions < z < 1.65 and < u; < 0.0005. The 
wave number is infinite at z = but decreases with an increase in z and increases with the increase in the 
angular frequency. The group velocity is greater than the phase velocity which indicates that the dispersion is 
anomalous. The value of refractive index is less than one which confirms this indication. 

The imaginary part of the determinant gives a dispersion relation quintic in k, i.e., Ai{z)k^ + A2{z,uj)k'^ + 
Arj{z,uj)k^ + A4{z,uj)k'^ + A5{z,uj)k + Aq{z,uj) — 0. This equation cannot be solved analytically for the exact 
solutions and is solved numerically using software Mathematica. The roots are approximated in arrays for each 
point of the two-dimensional mesh with equal step lengths 0.2 for z and lu. This equation gives one real root 
and two sets of complex conjugate roots. The real data values for the root gives a real interpolation function. 
This function represents the value of k shown graphically which is used to investigate the properties of medium. 

When we substitute the assumptions, the approximated root becomes infinite at z = due to which no wave 
exists there. We omit the value z = and our mesh reduces to 0.2 < z < 10, < < 10 for the interpolating 
function. The corresponding real root is shown in Figure 8. 

Figure 8 shows that the wave number is high near the horizon. It abruptly decreases and gradually increases 
as we go away from the horizon. The real wave number is shown in the figure. The wave number becomes 
negative in the region 0.845 < z < 10. Also, the phase and group velocities as well as refractive index are 
negative there which is the property of Veselago medium. The phase velocity is greater than the group velocity 
in the whole region admitting real waves. Also, the change in refractive index with respect to the angular 
frequency is greater than zero for the region 1.61 < z < 10, < uj < 10, thus it admits normal dispersion of 
waves [341 . The change in refractive index with respect to the angular frequency is less than zero for the region 
0.2 < z < 1.61, < w < 10 which shows anomalous dispersion of waves. 

6 Conclusion 

We derive the GRMHD equations for the Schwarzschild black hole magnetosphere in 3-1-1 formalism. We discuss 
one dimensional perturbations in perfect MHD condition. This has been explored for isothermal state of plasma. 
These equations are written in component form and then Fourier analyzed by using the assumption of plane 
waves. The determinant of the coefficients is solved for the following three cases in the magnetosphere: 

1. Non- Rotating Background (either non- magnetized or magnetized); 

2. Rotating Non-Magnetized Background; 

3. Rotating Magnetized Background. 

The non-rotating background is the pure Schwarzschild geometry outside the event horizon. The rotating 
background is the restricted Kerr geometry in the vicinity of the event horizon which admits a variable lapse 
function with negligible rotation. The dispersion relations and consequently, the wave numbers are obtained 
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Surface Plot of Wave Number 




Figure 6: The wave number decreases as the waves move away from the event horizon. The phase and group 
velocities are equal in most of the region. The dispersion is not normal. 
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Surface Plot of Wave Number 




Figure 7: The waves decreases in number as we move away from the event horizon of the black hole. The group 
velocity is greater than the phase velocity. The dispersion is anomalous in the region. 
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Surface Plot of Wave Number 




Figure 8: Some of the region possess the properties of Veselago medium. The phase velocity is greater than the 
group velocity in most of the region for which the dispersion is normal. 
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for these backgrounds. The properties of plasma are inferred on the basis of this number and the corresponding 
quantities are obtained in the form of graphs. 

We find that the gravity effect can be observed in the perturbations both in the non-rotating and rotating 
backgrounds. The values of k only show the gravitational effects on the smooth harmonic wave type pertur- 
bations. The rotation of the background adds this effect to the perturbed quantities. It can also be seen that 
the magnetic field perturbations are influenced in the case of rotating background whereas there is no effect of 
gravitation on the perturbed magnetic field in the non-rotating background. This implies that rotation affects 
the magnetic field. It is interesting to note that wave numbers decrease and the phase and group velocities 
increase as we go away from the horizon in each case. 

Recently, we have discussed the properties of the Schwarzschild magnetosphere for the cold plasma case 
[29j and some interesting consequences have been obtained about the wave analysis. In this paper, we have 
extended this analysis by adding the pressure effect (isothermal plasma). It turns out that pressure has brought 
changes in the case of non-rotating, non-magnetized background. In the case of cold plasma, normal dispersion 
of the waves was found for the non-rotating background whereas in the similar case of isothermal plasma, the 
dispersion of waves is not normal. In the case of rotating non-magnetized background one case of normally 
dispersed waves was found whereas in isothermal plasma with constant pressure, no such case is observed. This 
means that the background pressure ceases normal dispersion of waves. 

In the case of rotating magnetized background, there is one case of normal dispersion and two of anomalous 
dispersion. The case with normal dispersion of waves occurs when the plasma admits the properties of Veselago 
medium whereas in cold plasma there were two such cases where the plasma shows the properties of this medium. 
Thus the cold plasma equation of state provides more chances to the fluid to admit the properties of Veselago 
medium. 

It is worth interesting to note that the wave number becomes infinite and consequently waves vanish at the 
event horizon in all the cases. This corresponds to the well-known fact that no information can be extracted 
from a black hole. The negative phase velocity propagation regions are found in the cases of rotating background 
which are discussed by Mackay et al. in [30 according to which rotation of a black hole is required for negative 
phase velocity propagation. The same results were found in the case of cold plasma. Thus it can be deduced 
that the rotation brings negative phase velocity whether pressure is involved in the surrounding plasma or not. 
Mackay et al. |35j found that negative phase velocity propagation was neither observed in Schwarzschild-anti de 
Sitter spacetime, nor found outside the event horizon of a Schwarzschild black hole due to zero rotation. This 
is what we have obtained for non-rotating case. Ross et al. [36] also investigated the negative phase velocity 
propagation phenomenon for Kerr-Ncwmann and Kerr-Sen metrics particularly close to the outer event horizon 
when the magnitude of the charge is large. 

The complex solutions of the dispersion equation are not included because it is felt that these solutions 
will prove to be of little significance compared to those which have been discussed. The dispersion relations 
are obtained using 3-1-1 ADM formalism and contain the factor of acceleration (depend on lapse function and 
equals to —g) which makes them different from the usual MHD dispersion relations. We have investigated the 
waves propagating in a plasma influenced by the gravitational field. These internal gravity waves interrupt the 
MHD waves, therefore we have cases of dispersion in each hypersurface. It is mentioned here that some of the 
figures have patches missing which are due to the existence of complex values there. Mathematica cannot plot 
complex numbers with real numbers. These complex numbers are shown by the gaps in the figures. It would 
be interesting to investigate this analysis for the Kerr spacetime. Presently, this work is in progress. 
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Appendix A 

This Appendix includes the details to reach at the perturbed form of the GRMHD equations meant for the 
non-rotating background. The component form of these equations is also given. 
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The perturbed flow of isothermal plasma shall be characterized by its fluid density p, pressure p, velocity 
V and magnetic field B (as measured by the FIDO). The first order perturbations in the above mentioned 
quantities are denoted by Sp^ Sp, 5V and (5B. Consequently, the perturbed variables take the following form 

P = P° + Sp = p"+pp, 

P = p" + 6p — p'^ + pp, 

B = B" + (5B = B° + Bh, 

V = V° + (5V = V° + v, (Al) 

where p*^, p*^ , B*^ and V° are unperturbed quantities. The waves can propagate in 0-direction due to gravitation 
with respect to time t, hence the perturbed quantities depend on z and t. 

For the non-rotating background, the perturbed flow of fluid is only along z-axis. Thus the FIDOs measured 
four-velocity and magnetic field are along z-axis which can be expressed as V = u(z)ez, B = B{z)ey, and the 
Lorentz factor is 7 = ^^=5 • The perturbed quantities admit the following notations 

p= — ^p{t,z), p = — ^p{t,z), V = SY ^ v^it, z)e^, h = ^ ^ bz{t, z)e^. 
p p B 

Introducing perturbations given by Eg. ljAip in the GRMHD Eas. (|2.1ip - (|2.15p . we obtain 

= V X (av X B) + V X (aV X SB), (A2) 
V.{(5B) = 0, (A3) 
^feM + (.V.V)(^P + M + (P + P)7^V.| 

-a{p + p)yr.V\nu + a(p + p){W.v) + {Sp + 5p)V .{aV) + {Sp + Sp)j^Y.{aY.V)'V 
+2{p + p)j^{V.v){aV.V) ln7 + (p + p)j^{aV.VV).v + {p + p)j^V.{aV.W)v ^ 0, (A4) 
X 9 , ,4 1 ,ldv^ 1, , \d{SB),, 

+ {P + Ph^^^.3y' + {p + Ph^V,Vj,kV^V'' - -^{{aSB,),j - {a6B,)4B^ 

1 

-ip + p)lHvrV^ +v^VdykjV'' ~-/^{{Sp + Sp)V^ +2{p + p)-f^{Y.v)V^ + (p+p)w^}F„- 

--/^V^{{Sp + Sp)V^ +4{p + p)-f\V.v)V^ + {p + pyjVj^kV'', (A5) 

7' - ^^^^^ + - (P + P)7' V. ^ - 2 (p + p)7^ (V.v) ( V. V) In + 6(p + p)7' ( V.v) { V. ( V. V) V} 

a ot a at 

+ {p + p)74V.(v. V)V + 2(p + p)-f^V.(V.\/)v + 2{Sp + Sp^f^V-a + 2{p + p)7^(V.v) (V.a) 

+ (/5 + p)7^(V.v) + {p + p)j^v.a + -f^{Y .V){Sp + Sp) - i^iM _ (p + p)72(v.V) Inw 

a ot 

+2{Sp + Sp)j^V.{V.V)V + {Sp + Sp)j\\7.V) + 2{p + p)j^{V.v){V.V) + 2{p + p)-f^v .{V .\7)V 

-^-^[(V X B).(V X {aSB)) + (v x B).(V x {aBj) + (V x (5B).(V x (aB)) 
47ra 

+ (V X B).(V X ^) + (V X B).{^ X B)] = 0. (A6) 
The component form of Eqs. (jA2p - (|A6[) is given as follows 

a ot 

bz.z = 0, (A8) 

dp , _ _ , dp 2 / \dvz , s/-, 2 2\ 
P^+y-a{pp,z+PP,z) +P^+1 u{p + p) — + a{p + p){l + u )w^,2 

u' 

-{P ~ p){{oiup)' + aj'^u'^pu'} = a{p + p){l - 27^u^)(l-|- 7^-u^)— i;^, (A9) 



-72 {(dp + Sp) + 2(p + p)72(V.v)}a, - {Sp),, + —{{aB,),, - {aBj)4SB^ 
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1 dv 

{P + P)7^(l + l^u') i — 7^+ uv,^, } =~{pp + pph'^ia, +u{l+ j'u^)u'} 



a dt 

-l'^{p + p){u'{l + 7^"^)(1 + 47^^) + 2u'y^a^)}v^ ~ pp' ~ pp' , (AlO) 

nl d ^ 2^0^ 2 4 dVz I dp , ^ ~\ 2 r, 2 2 / n 

--^P + P^ -7nP+ -{p + Pn ~P--^ + [PP + PPn i^ua^ + 2-f u u + u \ 
a ot a at a ot a at 

+j^u{p,zp + p^zP + P.zP + P,zP) + {p + p)7^(l + 27^M^)u^_^ 

+ (P + P)7' |(37'mw' + a,)(l + 272^2) _ ^1 „^ = 0. (All) 

The conservation law of rest-mass [21] in three-dimensional hypersurface for isothermal state of plasma a{p + 
p)"fu = constant is used to obtain Eqs. (|A6|) and (IA9p . 



Appendix B 

When we consider non-magnetized isothermal plasma in rotating background, i.e., B = 0, the GRMHD 
Eqs. (|2.1ip and (|2.12p vanish and Eqs. (|2.13p - (|2.15p change into general relativistic hydrodynamical (GRHD) 
equations given as follows 

^^^tP) + (aV.V)(p + p) + (p + P)7'V.^ + (p + p)72v.(aV.V)V + (p + p)V.(aV) = 0, (Bl) 

{p + p){l^h+l'y^yo)] f-l: + V.v) = -(p + p)7'a, (B2) 



dt 



^'(-^ + ^-^^^P + P) + + P)7'V.(i| + V.V)V + 2(p + p)72v.a 

a ot act a at 

+(P + P)7^V.V = 0. (B3) 

We assume that the fluid is moving in xz-plane hence the fluid measured four-velocity by FIDO can be 
written as V = y(z)ex 4- u{z)ez for which the Lorentz factor becomes 7 = — \ . The following notations 
will be used for the perturbed quantities 

P= — =p{t,z), p= — ^p{t,z), V EE (5V = Z)ex -I- 2)62;. 

P P 
Introducing perturbations in Eqs. (|Bl[) - (|B3|) . it follows that 

^i^P±M. + (p+p)^2y^^ ^ a{p + p)j'V.{V.\7)^ + a(p +p)(V.v) 
+ {6p + 6p)\7.{aV) + {6p + Sp)-f^V .{aV .V)V + {aY.V){Sp + dp) 

+2a{p + p)j^{V.v){V.V)lnj + a{p + p)j'^{V.VV).v-a{p+p){v.\/\nu) = 0, (B4) 

1 dv^ 
a dt 

-^mSp + 6p) + 2{p + p)72(V.v)}a, - {5p),, - {p + p)^\v,V' + v'V,)Vk^,V'' 



{(p + p)7%- + {p + p)^^V,Vj]—^ + {p + p)l^v,,jV^ + {p + p)i^V^v,,kV'V'' 



-^'^{{5p + 5p)V^ +2{p + p)-f'^{W.^)V^ + {p + py}V^^J 



-7V,{((5p + dp)V^ + 4(p + p)YiY.w)V' + {p + pyjVj^kV", (B5) 

^2l d{dp + dp) ^ 2 ^ p)^4v.^ _ 2{p + p)74(V.v)(V.V) lnu + 6(p -I- p)7'(V.v){V.(V.V)V} 
a ot a Ot 

+ {p + p)j^V.{v.\/)V + 2(p + p)7'*{V.(V.V)v + v.(V.V)V} + 2{5p + 5p)-i^\.& 

+2(p + p)74(V.v)(V.a) + (p + p)72(V.v) + {p + pyi''w.Si + ^''{^ .V)[5p + 5p) - i^M 

a ot 

-{P + p)7^v.V) Inu + 2{Sp + (5p)7^V.(V.V)V + {Sp + 6p)j^{V.Y) 

+2(p-hp)7'*(V.v)(V.V) = 0. (B6) 
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The component form of Eas. (|B4|) - (|B6|) is 

dp , ^dp , ^„ , 2 /t/'9^^x (9Wz I . - , / , ^r 2t/ , /I , 2 2^ 



--(/5 - P)(au7p),2 + q;(p + p)7^w{(l + 27^1/^)^' + 2-f'^uVu'}vx 
7 

-a(p + p)|(l-27V)(l+72u^)^-27VV^t/'|i;, =0, (B7) 

1 CJU 1 CJU 

{p + ph\l + 7^2)--^ + {p + p)^\V-—^ + (p +p)72u(l + 7^2)1;,., + (p + p)7VF«,,, 

= -{pp + pp)-i^u{{\ + 7^2)1/' + 7\Fu'} - (p + p)7''u{(l + 472y2)uu' + WV'{\ + 7^ V"^)K 

-(p + p)7^[{(l + 2-i'^u'^){l + 27^2) - 7^2}]/' + 27^(1 + 27202)^^^^']^^^ (B8) 

= —{pp + pp)72{a2 + (1 + "f'^u'^)uu' + 7^^^^^'} - p'p - pp^z 
-{p + p)[l'^{u^V'{\ + 472^2) + 2V{az + (1 + 272u2)W}}]w^ 

-(p + p)7^[w'(l + l^u^){l + 472^2) + 2M72{a, + (1 + 2-i^u^)VV']]vz, (B9) 

7'-f + (7' - + -(P + P)7' + + 2(p + p)7^«l^«.,. + (P + P)7'(l + 27^u2)^;,,, 

Of OT a at a v ot dt/ 

+p[2p-f'^u{az + ^"^{yV + W)} + 72mp'] +p[2p72M{a2 + ^i^iyV' + W)} + 72-up'] 
+^'^upp^z + ^"^u pp + 7^Mpp.z + ^"^u'pp + 2(p + p)7'*{372wT^(uu' + FF') + uFa^ + uV}vx 

+ (p + p)7^ (67''w2(T/T/' + W) + -i'^{VV' + uu) + a;,(l + 27^712) - — + 27^^ I = 0. (BIO) 



u 



Appendix C 

The GRMHD equations for the rotating magnetized background with isothermal state of plasma remain the 
same as given by Eqs. p.lip - (j2.15p . Hence their perturbed form will remain the same as given in Appendix A 
(i.e., Eqs.dm-JMl)). 

In this case, fluid's four-velocity is the same as given in the previous section. The rotating magnetic field 
can be expressed in xz-plane, i.e., B = B[\{z)e^ + e^. We shall use the following notations for the perturbed 
magnetic field in addition to the notations given by Eq. (14.11) 

h=^ = b.j:{t, z)ex + bz{t, z)e^. 



Hence, the component form of Eqs. (IA2p - (IA6p can be written, after a tedious algebra, as follows 
dbx 

+ aubx,z = a'{vx - Xv^ + Vb^ - ub^) + a{vx.z - Xvz.z - X'vz + V'bz + Vbz.z - ub^), (CI) 

dh 

+ aubz.z = 0, (C2) 
hz,z = 0, (C3) 

dp dp , , 2 frr^^'^x dVz\ _ ^ f \c 2ta / 2 2n 1 

P'^ + + [P+ Pn \ + u—- > + aupp^z + aupp^z + ct{p + p){7 Vuv^^z + (1 + 7 " I'^z^z} 



dt dt t dt dt , 

--{p - p){au'-fp) z + a(p + p)72M{(l + 2j'^V'^)V' + 2^'^uVu'}vx 
7 

-a{p + p)\{l~2-f'^u^){l+-l'^u^)- -2-f^u^VV'\vz =0, (C4) 



u 

2n , 2,,2^ B^\ldu. r , . 4 Ai?2] IdVz 



(p + p)72(l + -f'V') + — ) — T^ + {{P + p)TuV 



An \ a dt 1 An \ a dt 
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B 



4tt 



XB^ 

47r 



B 
47r 



B' 
Aira 



bx{a'{l - u^) - auu'} + {pp + pp)j^u{{l + j'^V^)V' + -/^uVu'} 



{p + p)i'^u{{\ + AYV^)uu' + 4(1 + YV^)VV'] + 



B'^ua' 
Ana 



B^u. 



(p + p)7 [{(1 + 27V)(1 + 27^y^) - YV^}V' + 2Y{1 + 2Yu^)uVu'] + (Aa)' 

Ana 



vz = 0, (C5) 



[P + Ph\i + I'u') + ^1 -% + Up + ph'uv 



B' 



An j a dt 
X^B^' 



XB^ ] 1 dvx 



An 



,z + { {p + p)ruV + —) 



An I a dt 
XB^ 



+ ^{-{aXy + a'X - uX{uay}bx + (pp + pp)-f^[a, + u{{l + -f^u'^)u' + j^uVV'}] 
Ana 



{p + p)T{u^V'{l + AYV^) + 2V{a, + uu'{l + 2jV))} + 



XB^a'u 



Ana 



XB^u 
Ana 



{axy 



{p + p)tW{1 + 7 " )(1 + 47 ^i ) + 2uY{{l + 2Yu^)VV' + a,}} 

2 ^dp 2 1^ , \ T/N/ T/^^^ , \^"^ ^"^N 

^'am^^^ - 'Km + 4^^"^ - ~ + - ^) 

2 di) dv 

+ -ip + ph\V—^ + u—^) + 2{p + p)^^Vuv,,, + (p + p)7'(l + 2-/^u^)v,^, 
a at at 

+p^'^[2pu{a^+"f'^{VV' + uu')} + p'u + pu'] + p'-f'^[2pu{a^ + j'^{VV' + uu')} + p'u + pu'] 
B'^V 

+j'^upp^z + j'^upp.z - {aXyb^ 



Ana 



B^u. ,,,, B\uX-V). ^ 

-{aXybx + — ^ -{ab^.z + a'bx) 



Ana 



Ana 



(p + p)Y i &Tu^{VV' + uu') + Y{VV' + uu') + (1 + 2Yu^)az + 27^^ } + (aA)' 



B^X, 



(C6) 



Ana 



+Vx 



B 

{p + p)7'*{2MFa^ + Q-i'^uV{VV' + uu') + 2uV'} (qA)' 

Ana 



^ 0. 



(C7) 
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